Intrinsic spin-Hall accumulation in honeycomb lattices: Band structure effect 
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Local spin and charge densities on a two-dimensional honeycomb lattice are calculated by the Landauer- 
Keldysh formalism (LKF). Through the empirical tight-binding method, we show how the realistic band struc- 
ture can be brought into the LKF. Taking the Bi(l 1 1) surface, on which strong surface states and Rashba spin- 
orbit coupling are present [Phys. Rev. Lett. 93, 046403 (2004)], as a numeric example, we show typical intrinsic 
spin-Hall accumulation (ISHA) patterns thereon. The Fermi-energy-dependence of the spin and charge transport 
in two-terminal nanostructure samples is subsequently analyzed. By changing Ef, we show that the ISHA pat- 
tern is nearly isotropic (free-electron-like) only when Ef is close to the band bottom, and is sensitive/insensitive 
to Ef for the low/high bias regime with such Ef- With Ef far from the band bottom, band structure effects 
thus enter the ISHA patterns and the transport direction becomes significant. 

PACS numbers: 73.23.-b, 72.25.-b, 71.70.Ej 



CD 

S 

a 

B 

I 

o 
o 



> 

rn 
O 

o 



X 



In electron systems, the extrinsic spin-Hall effect (SHE), 
theoretically proposed long time ago,'-^ has been experimen- 
tally proven optically in semiconductor bulk structures^ and 
two-dimensional electron systems (2DESs)^ and even elec- 
trically in diffusive metallic conductorsi^ Recent achievement 
on observing the SHE in 2DESs at room temperatures* has 
further enlightened the possibility to manipulate spins via the 
SHE based on such a simple mechanism: transverse spin sep- 
aration by passing through longitudinal electric cuiTents. 

As for the intrinsic SHE, theoretically proposed much later 
than the extrinsic one^^ experimental evidence for its existence 
has been achieved only in two-dimensional hole systems^ but 
not in 2DESs. In particular, the local spin scanning in real- 
space for intrinsic SHE systems is difficult to carry out due 
to its limited size, and hence the required extra high reso- 
lution. Contrary to the extrinsic type, the intrinsic spin-Hall 
accumulation (SHA) pattern shows not only the out-of-plane 
component of spin accumulating antisymmetrically at the two 
lateral sides near the sample edges, but also oscillations due 
to the wave function modulation. Moreover, the SHA pattern 
may vary with, e.g., bias strength, spin-orbit coupling (SOC) 
strength, sample size and shape^^ 

In this paper we further investigate the crystal-structure- 
dependence, and hence the band structure effect, of the in- 
trinsic SHA (ISHA) pattern in 2DESs. The honeycomb lattice 
structure is particularly suitable for such investigation due to 
its nontrivial band structure and interesting geometry. In ad- 
dition, recent confirmation of the strong Rashba SOC" on Bi 
surfaces,'^ in particular the (111) case,'^ in which the pro- 
jected bilayer structure exactly forms the honeycomb lattice, 
provides a good numerical example to demonstrate these ef- 
fects. Based on the Landauer-Keldysh formalism (LKF),'^ 
the local spin densities (LSD) on four-terminal nanostructure 
samples made of the honeycomb lattice [see Fig. [TJa)]. are 
calculated. We will show that the ISHA pattern (i) exhibits 
isotropic/anisotropic spin transport behaviors when the Fermi 
level is near/away from the band bottom, (ii) shows dramatic 
difference between the left-right (zigzag) and the bottom-top 
(armchair) transport modes when the Fermi level lies in the 
band gap, and (iii) is extremely sensitive to the Fermi level in 



the low-bias regime. 

To apply the LKF, the first step is to construct the real space 
tight-binding-like Hamiltonian, which builds the underlying 
band structure and therefore is decisive for the transport prop- 
erties. For the honeycomb lattice and considering both ki- 
netic and Rashba hoppings up to the nearest neighbor only, 
the model Hamiltonian can be written as'^ 



itR{a X dj,),,] Ci, (1) 



with cj (ci) being the creation (annihilation) operator of the 
electron on site i. Here the first term is the on-site energy 
with parameter e^, which may also describe the local poten- 
tial and the disorder. In our potential-free case here, Si simply 
corresponds to the band energy offset in the language of tight- 
binding model (TBM). The second term in Eq. ^ contains 
the kinetic and Rashba hoppings with strengths to and tji, re- 
spectively. In the Rashba hopping term, dji is the unit vector 
pointing from site i to site j, and a = (cr^, a"^ ,a^) is the Pauli 
matrix vector. 
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FIG. 1: (Color online) (a) Schematic of the four-terminal setup with 
the conducting sample made of honeycomb lattice structure, (b) 
Band structures calculated by first principle (red dots) and by em- 
pirical TBM (black lines). 
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To obtain realistic values for the above parameters Si, to, 
and tji, we next perform the empirical TBM'^ of the Slater- 
Koster type,'^ based on Eq. ([TJ. Considering the two triangu- 
lar sublattices forming the honeycomb lattice and taking only 
single orbital pz on each site into account, it can be shown that 
the Hamiltonian matrix equivalent to Eq. ([U reads 

with the diagonal element given by Hn — Epl2, where Ep 
and I2 are the p-orbital energy and the 2x2 identity matrix, 
respectively, and the off-diagonal element given by 

„ U{1 + 2F) -^tn{l~F~V^G)\ 

^ -ztfl (1 - F + V3G) U{1 + 2F) J' 

(3) 

where U = l^Yppa + (l — ll) VppTr is the two-center inter- 
action integral involving pz atomic orbitals, and the compact 
functions are given by F = exp{—i^/3kya/2) cos{kxa/2) 
and G = exp{—i^/3kya/2)sm{kxa/2). Note that for 
an ideal (flat) two-dimensional honeycomb lattice such as 
graphene, only the tt bands contribute to U due to the van- 
ishing direction cosine 1^ = 0. However, later we will choose 
the Bi(lll) bilayer structure as a numerical example, in which 
both (T and tt bands contribute due to the nonvanishing l^- 
Still, here in both cases only the composite parameter U is 
to be extracted, instead of the individual Vppa and V^p^r ■ Ex- 
cept the Rashba hopping strength which we denoted con- 
sistently in both Eq. ([Ij and the empirical TBM Eqs. (|2]l and 
(|3]l, the other two parameters are related by —to — U and 
£i = Ep. Diagonalizing the 4x4 matrix of Eq. (|2]l gives the 
two pairs of the four energy dispersion curves. 

To extract reasonable parameters for Eq. ([1]), we consider 
the Bi(lll) bilayer structure, for its strong surface states,'^ 
making the electron transport 2DES-like, and its strong 
Rashba SOC,'^ making the ISHE thereon promising. Con- 
sidering nearest-neighbor hopping only (and hence neglect- 
ing the interbilayer hopping), the projected two-dimensional 
lattice structure is exactly of honeycomb type. We therefore 
fit the energy dispersion curves obtained from diagonalizing 
Eq. with the surface band structure from the first principle 
calculation,'^ as shown in Fig. [Tfb). Both directions in this 
plot are along FA/. Due to the simple model we have taken, 
the fitting gives good agreement only near the F point. Cor- 
respondingly, the band parameters extracted from such fitting 
areio = 1.6302 eV, tR = 0.1853 eV, ande, = 4.8324 eV. 

The lead-sample setup is sketched in Fig. [TJa). where the 
sample made of the honeycomb lattice is of nearly square 
shape with four terminals (left, right, bottom, and top), each of 
which can be contacted by a semi-infinite normal metal lead. 
Throughout the rest of the calculations, we will set the x and 
y axes along the zigzag and armchair directions, respectively, 
and fix the origin of {x,y) = (0, 0) at the center of the sample, 
as indicated in Fig. [Tfa)- 

The LKF, namely, the nonequilibrium Keldysh Green's 
function formalism'^ applied on the Landauer multiterminal 
setups, has been summarized in some detail in Ref. |9t Its 
central spirit is to solve the kinetic equation for the Keldysh 
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FIG. 2: (Color online) Left panel: the band structure along VM di- 
rection calculated by the empirical TBM. Right panel: the total den- 
sity of states of the sample with size about 4.1 x 4.2 nm'^ calculated 
by the LKF. The dashed lines in both panels are guided for the eyes 
and the dotted lines in the left panel indicate the energy range shown 
inFig.db). 



Green's function matrix G"^ {E), written in the real-space rep- 
resentation at energy E. Physical quantities such as local 
charge density (LCD), LSD, or even the charge and spin cur- 
rent densities, can then be extracted from {E). In systems 
free of phase-breaking interactions (such as electron-electron 
or electron-phonon interactions), the self-energy modifying 
the carrier life time inside the sample is contributed only from 
the contact leads through the nearest neighbor hopping with 
strength set equal to to [see the bold (red) connection lines in 
Fig. [Ila)]> and the Keldysh Green's function can be solved 
exactly.'** Explicitly, the self-energy can be expressed as the 
product of and the retarded surface Green's function of the 
attached semi-infinite leads/.iSiL'^ In our analysis, we will con- 
centrate on z-component of the LSD given by^ 

t l-EF+eVa/2 

(^M>r. = Tr [a'^G< (F; r„ r,)] dE, (4) 

■^JEF-eVa/2 

where = {h/2) is the /i component spin operator with 
^ — x,y,z and is the position vector of the ith lattice site. 
In the right-hand side, the 2x2 matrix {E; r^, r^) is the ith 
diagonal submatrix element of the whole {E), Ep is the 
Fermi energy to be tuned in the later analyses, e is the electron 
charge, and Vb is the applied potential difference between the 
negatively and positively biased leads. The LCD, which will 
be shown to modulate the LSD in the later analysis, is given 
by 

(e7Ve)r, =e/ TrG< (F;r„r,)dF, (5) 

JEF-eVo/2 

where is the electron number operator Note that here 
the band bottom F^, which is set to the lowest energy of the 
full band structure calculated from the previously introduced 
TBM, does not explicitly enter the expressions (|4|i and (|5]) but 
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FIG. 3: (Color online) Local spin density on a 2.3 x 2.2 nm sample 
with low bias with the Fermi level Ef set at (a) 0.08 eV above the 
band bottom, (b) the middle of the lower band, (c) the top edge of 
the lower band, and (d) the middle of the energy gap. The red/green 
(dark/light) dots represent positive/negative (Sz) with the dot size 
proportional to the magnitude. The maximum values of (Sz) are of 
the order of 10"'^ {h/2) in (a), and lO""* (ft/2) for (b)-(d). 



will play a decisive role in dealing with the self-energy due to 
the lead coupling. It may become even more crucial when cal- 
culating physical quantities to which equilibrium states also 
contribute, such as the bond spin current density.^ 

Before performing the LSD and LCD calculations, one 
last step is to report the convincing correspondence between 
the band structure of the infinitely extending honeycomb lat- 
tice, calculated by the empirical TBM, and the total density 
of states (DOS) of the finite-size sample, calculated by the 
LKF. The latter is given by summing the spectral function 
for each site pr (E) = -i Tr {im [G^ (r^, rf,E)] } /tt, 
where G^ (ri,ri;£') is the spin-resolved ith diagonal sub- 
matrix element of the retarded Green's function matrix. As 
shown in Fig. |2] the main features, including the band top, 
band gap, and band bottom, are nicely correspondent with 
each other. Note that the drop of the total DOS at the band gap 
region cannot be perfectly step-function-like since the sample 
considered in the LKF is not infinitely large. 

Having constructed the correspondence between the empir- 



ical TBM and the LKF, and extracted reasonable parameters 
by fitting to the first principles calculation for the Bi(ll 1) sur- 
face, we are now ready to present the local spin, and later also 
charge, densities on the conducting sample in the honeycomb 
lattice. Two bias regimes will be distinguished: low and high, 
standing for bias voltages of eVo = 0.002 and 0.2 eV, respec- 
tively, in the rest of the paper Although there are totally four 
terminals free to contact the electrodes, we will consider only 
two-lead cases with head-to-tail orientation, either parallel to 
the X or y axes. 

We first present the LSD for a sample of totally 248 lattice 
sites (sample area about 2.3x2.2 nm^) in the low bias regime. 
To examine the direction dependence of the spin transport, we 
set E'i? at some representative positions. As shown in Fig. |3] 
we gradually raise Ep from the band bottom (Ep — 0, i.e., 
about 0.083 eV above the bottom of the band, consistent with 
the first principle calculation) to the middle of the band gap 
(about Ef = 4.82 eV). In each panel of Fig. |3] the {S^) 
distribution is antisymmetric about the bias axis, and there- 
fore exhibits the main feature of the intrinsic SHE. Further- 
more, as can be clearly seen that, only when Ep is set at the 
band bottom [Fig. |3la)], the pattern becomes nearly isotropic, 
i.e., the spin transport does not show direction dependence 
and thus behaves as free electrons. With the increase of Ep, 
Figs. [3jb) and |3jc) show distinct LSD distributions for the 
two different transport modes. Interestingly, when Ep is set 
in the middle of the gap [Fig. [3jd)], a dramatic difference 
between the two transport modes emerges. In the left-to-right 
case spin accumulation occurs at the lateral edges, while in the 
bottom-to-top case there is almost nothing at the lateral edges 
but some spots induced near the leads. This can be understood 
by observing that at this Fermi energy the transport inside the 
sample is supported by the edge states, which are contributed 
mostly from the zigzag edges. Thus in the bottom-to-top ge- 
ometry the whole sample behaves nearly insulating: neither 
charge nor spin can pass through the sample. 

Next we take further looks at the spin and charge transport 
near Ep = 0. We consider a sample with 802 sites (about 
4.1 X 4.2 nm^) with two leads in the left-to-right orienta- 
tion, and finely raise the Fermi level from Ep — eV to 
Ep — 0.08 eV. Both the low and high bias regimes will 
be analyzed. As shown in the first column of Fig. ID the 
ISHA pattern in the low-bias regime is extremely sensitive 
to the Fermi level. Compared to the LCD in the second col- 
umn of Fig. m one can see that such sensitivity stems from 
the wave function modulation. Whereas the nonequilibrium 
transport is contributed from the states between Ep — eVb/2 
and Ep + eVb/2, the low bias regime with small voltages is 
mainly described by the Fermi energy state, and the electron 
transport thus behaves quantum mechanically. It turns out that 
the Fermi level, determining the length of the wave vector k, 
influences the formation of the electron wave, and hence in 
turn the ISHA pattern. It is interesting to note that a strong ac- 
cumulation of electrons does not necessarily lead to a promi- 
nent accumulation of spins. Conversely, where there are no 
electrons, there must be no spins. Put in another way, the lo- 
cal spin and charge density patterns must be, to some extent, 
consistent with each other. 
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FIG. 4: (Color online) LSD (the first and third columns) and LCD 
(the second and fourth columns) on a 4.1 x 4.2 nm^ sample. Each 
row corresponds to a gradually changing Ef near the band bottom. 
The left/right two columns are for low/high bias regimes. For LSD, 
the maximum value of (5*2) in each panel in the low and high bias 
regimes is of the order of 10"^ (h/2) and 10"^ (/t/2), respectively. 
For LCD, the maximum values of (A^e) in each panel in the low and 
high bias regimes are of the order of 10"* and 10""^, respectively. 



For the high bias case, both the local spin and charge den- 
sities change moderately with the increase of Ep (the two 
right columns in Fig. HJi. Contrary to the low bias regime, the 
involved states participating electron transport cover a much 
wider range of energy. Summation of these states with dif- 
ferent wave lengths eventually gives a waveless charge distri- 
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bution, as compared to the low-biased patterns. The electron 
transport behavior is therefore far from the standard quantum- 
mechanical description. In this case the nonequilibrium spin 
accumulation is no longer affected by the wave function mod- 
ulation, and the ISHA pattern is robust against the change of 
Ep. 

Before closing, it is worthy to remark here that in the pio- 
neering formulation of Ref. ^ and also the recent application 
of Ref. [1^ on the triangular lattices, the crystal structure in- 
formation is lost since the Fermi energy is chosen close to 
the band bottom: Ep — Ei, + 0.2to. Consequently, the spin 
accumulation property remains free-electron-Iike and exhibits 
rotational invariance (except for a coexistence of the Dressel- 
haus term, giving rise to anisotropic dispersion). It is only 
when the Fermi level is far from the band bottom, where the 
corresponding wave vectors are short, that the band structure 
(or the crystal structure) effect emerges. 

In conclusion, taking the honeycomb lattice as a particu- 
lar case, we have pointed out the crucial role Ep plays in 
the ISHE due to band structure effects. Recent observation 
of the strong surface state and giant Rashba SOC on Bi(lll) 
surface'-' has attracted much attention and is especially suit- 
able as a numeric example in our investigation. The possibil- 
ity of observing the quantum SHE on the Bi(ll 1) surface and 
its multilayer thin film''' has made Bi(l 11) even more promis- 
ing. Here we have reported another positive viewpoint of its 
potential of observing the ISHE thereon, provided that Ep 
and the transport direction are important. Moreover, we have 
combined the LKF with the first principle band calculation 
through the empirical TBM, allowing one to extract realistic 
band parameters for the LSD calculation. 
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